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I. INTRODUCTION 



The well-known complexities of 4D field theory have of- 
ten forced theorists to test models in lower-dimensional 
spaces. In general, the foundations of such models have 
been obtained only by projection from the physical di- 
mension which, of course, cannot shed light on the sub- 
tleties inherent to a particular dimension. In this vein — 
as it was pointed out by Binegar [1] with good reason — 
an independent development of the theories in their na- 
tive dimension is required since a theorist is not supposed 
to be omniscient. 

In this sense, three-dimensional theories deserve a spe- 
cial attention due to its closeness to reality. Fortunately, 
planar physics has undergone a remarkable development 
in the last few decades. A host of new experimental re- 
sults coming mainly from condensed matter physics and 
the accompanying rapid convergence of theoretical ideas 
have brought to the subject a new coherence and have 
also raised new interests. Among the so many and in- 
teresting planar models that have been investigated, it is 
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worth mentioning the graphene [2]. This genuinely pla- 
nar carbon system seems likely to be a good framework 
for the verification of ideas and methods developed in 
quantum (gauge-) field theories. Consequently, we hope- 
fully expect that the techniques of QED 3 when applied 
to this low-dimensional condensed-matter model lead to 
new and relevant results Hil. As far as gravity is con- 
cerned, the reason for doing research on planar gravity 
is quite amazing: (2+l)-dimensional gravity has a di- 
rect physical relevance to modeling phenomena that are 
actually confined to lower dimensionality. In fact, gravi- 
tational physics in the presence of straight cosmic strings 
(infinitely long, perpendicular to a plane) is adequately 
described by three-dimensional gravity |5|. We remark 
that the causality puzzles raised by 'Gott time machines' 
were solved with the help of this lower dimensional model 
0. 

On the other hand, three-dimensional electromagnetic 
(gravitational) models enlarged by a Chern-Simons term 
have been the object of much attention. In the vector 
case this term gives mass to the photon in a gauge- 
invariant way; while, for planar gravity the Chern-Simons 
term is responsible for the presence of a propagating 
parity-breaking massive spin-2 mode in the spectrum of 
the model Q. 

Recently, both higher-order electromagnetic and gravi- 
tational models have enjoyed a revival of interest. Indeed, 
the fourth-dimensional theory of quantum electrodynam- 
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ics proposed by Lee and Wick (LW) with the purpose 
of understanding the finite electromagnetic mass split- 
ting of mesons, prior to QCD was established [1, H|, has 
been rather explored as a kind of toy model for the more 
complex dynamics of the LW Standard Model, i.e., the 
model in which Grinstein, O'Connel and Wise, building 
on the pioneering work of LW, introduced non-Abelian 
LW gauge theories |10h19| |; whereas, just about three 
years ago, Bergshoeff, Hohm, and Townsend (BHT) [20I — 
|28I | proposed a particular higher-derivative extension of 
the Einstein-Hilbert action in three spacetime dimensions 
whose linearized version is a rare example of a fourth- 
order system that is not pestered by ghosts [29]. Besides, 
a canonical analysis of the quadratic curvature part of the 
BHT system done by Deser [29] establishes its weak field 
limit as both ghost-free and power-counting UV finite, 
thus violating standard folklore in the extreme. 

The preceding considerations naturally suggest that in- 
vestigations into general 3D higher-order electromagnetic 
(gravitational) models with a Chern-Simons term, are 
welcome. The introduction of higher-derivatives, nev- 
ertheless, could in principle jeopardize the unitarity of 
the models. It would thus be very convenient, in the 
spirit of paragraph one, to devise an easily handling pro- 
cedure, specific to planar models, which allowed, on phys- 
ical grounds, a constructive and meaningful discussion 
of the unitarity of generic 3D electromagnetic (gravita- 
tional) models, in an uncomplicated way. 

In this paper, the aforementioned procedure is con- 
structed by means of a new basis of spin operators, spe- 
cific to 3D models, which allows a Lagrangian decompo- 
sition into spin components. 

The ideas underlying our theoretical framework are 
described in Section HH We start off by building up a 
new class of spin projectors for 3D models and then 
discuss how to obtain the propagator for these models 
via the mentioned operators. The procedure for prob- 
ing the unitarity of the 3D models is constructed after- 
wards. Two important and timely higher-derivative sys- 
tems enlarged by both Chern-Simons and higher deriva- 
tive Chern-Simons terms are employed in Section IIIII to 
illustrate the level of generality and quickness of the 
method. In Section IV it is shown that the expressions 
"closure" and "completeness" cannot be used interchange- 
ably, as far as the projection operators are concerned, as 
it sometimes done. We also comment in this section on 
the wrong idea that Chern-Simons terms can be utilized 
for curing the nonunitarity of higher-derivative models. 
Some of the more technical results are gathered in the 
Appendix |A] 

In our conventions the Greek letters denote spacetime 
indices, the metric signature is (+1, -1, -1), £012 = +L 
where e^p is the Levi-Civita symbol, and h = c = 1 



II. PRESCRIPTION FOR PROBING THE 
UNITARITY OF 3D MODELS 

In the analysis of quantum aspects of any field theory, 
considerable interest is devoted to the description of the 
particle spectrum and the relativistic quantum proper- 
ties of scattering processes of the theory under investiga- 
tion. Some of these issues may be understood by means 
of the analysis of the propagator of the theory, which is 
obtained by the inversion of the wave operator. Accord- 
ingly, it is of great and fundamental importance to per- 
form this inversion judiciously. We shall begin by looking 
for a suitable basis for the linear operators acting on the 
fields of the model. Using this basis, a generic expression 
for the propagator will then be constructed. Finally, a 
procedure for analyzing the unitarity of the 3D models, 
based on the preceding ingredients, will be worked out. 

A. A new set of spin projection operators for 3D 
models 

We start off by searching for a basis for the vectorial 
space of the wave operators. The vector space where 
these operators act is formed by finite-dimensional rep- 
resentations of the Lorentz group. In 4D, for instance, 
it is always possible to decompose these vector spaces in 
a direct sum of subspaces with well defined spin since 
locally the Lorentz group can be regarded as the ten- 
sor product SU (2) ® SU (2). Besides, the only mapping 
operators that can be built among these projectors are 
those associated with the same spin. In fact, the ex- 
istence of mapping operators implies in a bijection be- 
tween spaces which can be achieved if and only if they 
have the same dimension. However, the construction of 
operators that map a spin space into a subspace with a 
different spin can only be realized by decomposing the 
larger spin-space into a direct sum of subspaces defined 
by preferential vectors. The explicit construction of the 
spin projectors for fields of arbitrary rank can be made 
by appealing to the tensor product of projectors of vec- 
tor fields according to the rules of group representations. 
The spin projectors, in 4D, that decompose the vector 
fields can be explicitly constructed using the Minkowski 
metric and partial derivatives, as follows 

P(0) ia/ = Ultv = ^-, U = d^, (la) 

P{^) fi u = e ^=ViMu-u) fJll/ . (lb) 

A careful analysis of the preceding equations allows to 
conclude that the spin projectors and mapping operators 
of spin subspaces with the same dimension should also be 
built solely with the metric and partial derivatives, which 
leads us to the well known Barnes-Rivers operators [30- 
[32j. It is worth noticing that if extra vectors are used 
in the construction of the models, such as a background 
vector in Lorentz violating models, operators with well 
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defined spin will be insufficient to form a basis for the 
wave operators [33[. 

On the other hand, the issue of the attainment of the 
wave operator and, subsequently, of the propagator, for 
3D models, needs to be dealt carefully. Why is this so? 
Because now we have both parity-conserving (PC) and 
parity- violating (PV) models (34|. Since the PC systems 
are defined by Lagrangians involving only the metric and 
partial derivatives, the appropriate basis for expanding 
the wave operator is, of course, that made up of the 
usual 3D Barnes-Rivers operators. The 'mark' of the 
PV models, i.e, the characteristic feature that enables us 
to recognize them, is, in turn, the presence of the Levi- 
Civita tensor, which allows us to define another vector 
linear operator, 



(2) 



Using this operator, we can enlarge the usual operator 
basis, in order to obtain a complete set of linear 

operators {9,uj,S} [35j |. It is worth noticing, however, 
that 9 is no longer a spin operator since in the massive 
case the 3D spin corresponds to unitary representations 
of SO (2), that are one-dimensional. In fact, the oper- 
ators 9 and uj divide the three-dimensional space into a 
direct sum of two subspaces with dimensions 2 and 1 , in 
this order, which implies that 9 does not project into a 
spin subspace. Consequently, it is impossible to put a 
transparent and accurate physical interpretation on the 
excitation modes related to the PV models if they are 
expressed in terms of the basis {9,to,S}. A successful 
way of dealing with this problem would be to opt for a 
basis associated with the spin of the particles in 3D. We 
discuss in the sequel how this basis can be constructed. 

We begin by recalling that in quantum field theory par- 
ticles are identified as unitary irreducible representations 
(irreps) of the Poincare group. This identification pro- 
vides two quantum numbers for particles: mass and spin. 
The spin is characterised by the unitary representations 
of the little group of the representative momentum of the 
particle, namely, the subgroup of the Lorentz group that 
leaves the representative momentum unchanged. 

In 4D, the task of identifying the spin operators is eas- 
ier than in 3D. The reason why this is so comes from the 
fact that in 4D the spin of the massive particles is given 
by unitary irreps of the group SO (3). Such representa- 
tions are associated with the representations of the group 
SU (2), which is the covering group of SO (3). Further- 
more, the fundamental representation of SU (2) is equiv- 
alent to its complex conjugate. Therefore all the rep- 
resentations of SU (2) are real and, consequently, repre- 
sentations of SO (3) are univocally related to the SU (2) 
representations. This implies that if the wave operator 
is decomposed into operators that projects into well de- 
fined irreps of SO (3), they are automatically identified 
as operators with well defined spin. 

In 3D, on the other hand, we have a distinct situation, 
since the unitary representations of SO (2) are associated 



with U (1). The fundamental representation of U (1) is 
not equivalent to its complex conjugate. Since all repre- 
sentations of SO (2) are real they are not directly related 
to representations of U (1), but rather they should be 
identified with the direct sum of the fundamental repre- 
sentation and its complex conjugate. However, all the 
irreps of SO (2) are two-dimensional and can be associ- 
ated with representations of U (1) by the complexification 
of the fields. Consider, for instance, the vectorial repre- 
sentation of 50(2). A transformation of SO (2) acting 
on a vector A = (Ai,A 2 ), yields 



cos# — sin6* 
sin(9 cos0 



A 2 



(3) 



The normalized eigenvectors of this transformation are 



Ai = 



V2 v« 



A 2 = 



V2\ 



(4) 



So, we can define a basis for the 3D Minkowski space, 
with the characteristic that each of its vectors spans one 
and only one ID subspace that is an eigenspace of the 
U (1) transformations, i.e., 



e(0) L 



e(+l)„ = (ei) u = 



e(-l) =(e 2 ) 




(5) 
(6) 
(7) 



where e (0) is time-like, whereas e\ and e 2 are space-like 
vectors. 

Thence, under a suitable unitary change of variables, a 
real vector field transforms like 



Ax 

-4-2 



A x 
A 2 



J_ / Ai + iA 2 
^2\Ax- iA 2 



(8) 



In this way the transformation of the vector A 1 under the 
rotation §3§ is given by 











Ax 
A 2 



(9) 



We conclude, therefore, that it is possible to make the 
identification of the vectorial representation of 50(2), 
A, as the direct sum of the U (1) fundamental represen- 
tation, □ , and its complex conjugation, □*, i.e., 



A - □©□* 



(10) 
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The 0-operator is the identity in the space of the repre- 
sentation A. We may split this space in the direct sum of 
one-dimensional subspaces. Keeping in mind the 3D spin 
representations, spin projection operators may be asso- 
ciated to the basic complex vectors e\ and e%, as follows 



T>"' 



(11) 
(12) 



Here, p is the projection operator associated with □, 
while a is related to □*. Note that p and a are related 
by complex conjugation, p* = <r; in addition, they are 
Hermitian and non-symmetric: 



Pp V = (ei) (ei)* = (e 2 )„ {e 2 )\ 



(13) 



This completes the identification of the spin modes for 
vector fields. It is also important to express the Chern- 
Simons operator @ in terms of the spin projection opera- 
tors. In order to do this, we note that £ A " yp e^e^-^= = —i, 

since ej., e 2 and are normalized, which allows us to 
write 



Sn 



k 2 (p^ 



(14) 



One could wonder about the possibility of building 
mapping operators between the subspaces defined by k 
and e\ or k and e 2 . In reality, these mapping operators 
are unnecessary since they would explicitly depend on e\ 
and e 2 . Actually, in the case of Lorentz preserving mod- 
els, the wave operator is constructed using solely 77's, d's 
and e's. 

Before going on, it is important to discuss the meaning 
of parity as far as the 3D operators we have just analyzed 
are concerned. In 3D, the representation of the parity 
operator in the Minkowski vector space is given by 



P 




(15) 



As a result, we get from flnj and ([12} that PpP^ 1 = a 
and PaP~ l — p, which clearly shows that, unlike what 
occurs in 4D, we cannot assign a definite parity for a 
given spin. We remark that 



PloP- 



P6P- 



(16) 



In the special case of the PC models, the aforementioned 
relations allow us to conclude that the massive excitation 
modes for nontrivial spins are nothing but spin doublets. 

After this little digression, let us build up the spin 
projection operator for rank-2 tensors. For these tensors, 
we have 



(le-ieo)® (ie 

(3 x 0® 2 x 1 



-1©Q) = 

2 x -10 20 -2), (17) 



where the underlined numbers denote the spin, and the 
remaining ones are related to the spin-multiplicity. 

Now, a general rank-2 tensor, T^ u , may be written as 
product of two vectors, say and B u . Therefore, 



PflU AnP>y 



(18) 



Since a generic vector can be split in its spin components 
A^ D (1 ffi —1 ffi 0), by means of the spin projection op- 
erators, p, a and w, namely, 



A, 



(Pup + a up + w w ) A p , 



(19) 



a generic rank-two tensor may also be decomposed in its 
spin components as follows 

Ppv (p ppPva ~t~ Ppp^va ffi Ppp^va ~t~ @ ppPvo ffi O pp®va 

+ evpuw + uJ^ppucr + u w o vcr ffi uj^uj^Tf" 7 . (20) 

Note that p, a, and uj are associated with spin +1, -1, 
and 0, respectively, which implies that pp, puj, up, pa, up, 
ww, auj, uja, 00 (with the indices omitted) are associated 
with spin +2, +1, +1, 0, 0, 0, -1, -1, and -2, in this order. 

In the case of the graviton field h^, which is a sym- 
metric rank-2 tensor, the symmetrization of the operators 
above yields the following spin projection operators 



P M (+2)„ = PppP 



P hH {~ 2 )pu;pa = °ppOvo 



(21a) 
(21b) 



P (+l) M „;p CT = 2 (Pm"*" + P»P U ^ + Ppo^vp + Pua^pp) , 

(21c) 

P {- 1 )p,u;pa- = 2 i a W U ™ + OvpUpo ffi (Tp,aU V p ffi <T ua U}p p ) , 



(21d) 

(21e) 



P 22 {°)p lV ;p < T - 2 (PvP a v<? + PvpVpo + Ppo^vp + PvoOpp) ■ 

(21f) 

The preceding operators, of course, are Hermitian. In 
addition, the projection operators associated with non- 
trivial spins are complex, whereas those related to spin-0 
are real because non-trivial spins are non-trivial represen- 
tations of U (1), that are complex. For a real Lagrangian, 
the complex structures (|21a[) - (|21dl) alone cannot appear 
in the wave operator decomposition in terms of the spin 
projection operators (this decomposition will be clarified 
later). We can ensure, however, due to the Lorentz invari- 
ance of the model, that projectors of the irreps of SO (2) 
will be present in the wave operator. Such operators, 
usually known as Barnes-Rivers operators, are written in 
terms of 9 and u). Using the identity, 9 = p ffi cr, we may 
split the real Barnes-Rivers operators into spin projection 
operators. Since p* — cr, the wave operator is obviously 
real. 

Consider, for instance, the projector associated with 
a symmetric and traceless rank-2 tensor. This operator 
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projects into a non-trivial and irrep of SO(2) and, there- 
fore, it is two-dimensional, and can be written as 

( 2 ) All y;p C r = ^ + ® p.c,Q V p) - ~0 p,y9 pa . (22) 

Now, taking into account that 9 = p + <r, we obtain two 
projectors in terms of p and a, one for each degree of 
freedom of spin, i.e., 

pHh ( 2 )^-pa = PwPvo + ewcw, (23) 

clearly showing that the Barnes- Rivers operator P hh (2), 
is nothing but a sum of spin +2 and spin -2 operators. 
This process of decomposition can be repeated for all 
operators needed to exhaust all the possibilities of con- 
traction of the fields present in the free Lagrangian. With 
this decomposition, the gravitational Chern-Simons op- 
erator 

Sfj,v;pa — 9ppS V( j -\- Qpa^up ~\~ QupSp jC! ~t~ 6 ^uo& "fipi (24) 

can be expressed as 

S„ = ~4zVP (P (+2)^ - P (-2)„) . (25) 

The other relations among the operators are listed in the 
Appendix [3] 

B. The propagator 

We are now ready to find the propagator and present 
afterwards the algorithm for probing the unitarity of 3D 
electromagnetic (gravitational) models. Consider, in this 
vein, a 3D Lagrangian C which is a function either of a 
vector field A a or of a symmetric rank- 2 field, h a (,. In 
order to compute the propagator for the model, we need 
beforehand the quadratic part of £, i.e., 

( £ )2 = 2 ^VaOaptpp, (26) 

a,P 

where a, j3 represent vectorial or tensorial indices, C Q( g 
is a local differential operator (the wave operator) and 
ip a encompasses the fundamental quantum fields of the 
model. For gravity models, for instance, this is accom- 
plished by means of the weak field approximation of the 
metric, i.e., = f]^ v + h^. 

Using the identities of AppendixjAj we then expand the 
wave operator in the basis of the spin operators, namely, 

Here, a (J)^ are the coefficients of the expansion of the 
wave operators. The diagonal operators, P^ v (J), are op- 
erators that project the field ip into its spin J. Whereas 



the off-diagonal operators (i ^ j) implement mappings 
into the corresponding spin doublet subspace. The re- 
sulting spin operators do obey the orthonormal multi- 
plicative rules and the decomposition of unity, i.e., 

Y^Piji^PkiiJ)^ = 5 jk 5 IJ P u (I) ail (28) 

P 

J2 P ™( J )*P=t a 0- (29) 

i,J 

This converts the task of inverting the wave operator (|27p 
into an straightforward algebraic exercise. Indeed, all we 
have to do is to invert the matrix of coefficients a ( J)ij- 
Nevertheless, a may be degenerate due to the gauge 
symmetries of the model since the physical sources actu- 
ally may satisfy some constraints. These consistently ap- 
pear in order to inhibit the propagation of non-physical 
modes. The explicit expressions for these constraints are 
given in terms of the left null-eigenvectors vj 1 "'^ of the 
degenerate coefficient matrices 

J2v} L ' n) (J)P kj (J) a0 Sp=O. (30) 

P 

Nonetheless, since the propagator is saturated with the 
physical sources, the correct procedure for the attainment 
of the propagator is to invert any largest non-degenerate 
(for general values of momenta k) sub-matrix of a (J)ij- 
Accordingly, in order to obtain the propagator, it suffices, 
in practice: (i) to delete rows and columns of a (J) t j ac- 
cording to the number of gauge symmetries, which gives 
rise to a matrix that we shall call A (J)y, and (ii) to in- 
vert A and subsequently saturate this matrix with 
physical sources. As a result, the saturated propagator 
(II) assumes the form 

Tl = iJ2S*aMJ)v P ii( J )<*f> S 'p- ( 31 ) 

It must be emphasized that since the physical sources 
satisfy the constraints (|30[) . the propagator is gauge in- 
dependent. That is a great virtue of our method in com- 
parison with the methods that do not use orthonormal 
projection operators. Indeed, in our procedure no gauge 
fixing is required. 

C. The prescription 

For the sake of simplicity, we shall divide our discussion 
about the unitarity of the 3D models into two parts: one 
of them related to the massive poles, the other concerning 
the massless ones. 

• massive poles 
To ensure that there are neither ghosts nor 
tachyons in the propagation mode of a given 3D 
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model, we must require that at each simple pole of 
the propagator (k 2 — m 2 ), 



3mRes(n| fe 2 =m 2) > 0, and m 2 > 0. (32) 

In the light of (13"TT) . we come to the conclu- 
sion that the condition for the absence of ghosts 
for each spin and for arbitrary sources is di- 
rectly related to the positivity of the matrices 

(E A (J, m 2 )y Pa {J)) ag , where A (j, m 2 )^ = 
Res A (J,m 2 )'. 



k 2 —m 2 



is the matrix A (J) ^ with 

the pole extracted. Furthermore, it can be shown 
that these matrices have only one non-vanishing 
eigenvalue at the pole, which is equal to the trace 
of A^ 1 (J,?n 2 )L 2 _ 2 . Besides, the own operators 
Pij (J) contribute only with a sign (— 1)^, when- 
ever calculated at the pole, where N is the sum of 
the number of p's and it's in each term of the projec- 
tor. This sign (— 1)^ can be understood by noting 
that in rest frame of the particle, p and a contribute 
with a minus sign, whereas lo contributes with a 
positive sign. In summary, we may say that the 
conditions for absence of ghosts and tachyons are 
such that for each massive single pole: (i) m 2 > 0, 
and (ii) (-l^trA- 1 (J,m 2 ) > 0. 

• massless poles 

The massless modes have some subtleties which re- 
quires an extra care. Indeed, at first sight it seems 
that the basis of operators is not well-defined for 
light-like momenta. However, the physical sources 
constraints, that have its origin in the gauge sym- 
metries of the model, allow to find a well-defined 
expression for the saturated propagator, even for 
light-like momenta. These physical constraints (|50|) 
take the form k^S^ = for the electromagnetic 
models and k^S^ = for the gravitational ones. 
Consequently, a convenient way of avoiding ghosts 
in the massless modes is to rewrite the inverse of the 
wave operator in terms of the following structures 



k^ k v 
~k?~ 



(33) 



This task can be greatly facilitated by appealing to 
the relations of the Appendix [A] in addition, the 
sources must be expanded in a suitable momentum 
basis, 



= cik,j, + c 2 g M + c 3 e M , (34) 

S^u = c 1 k il k l ,+c 2 {k l _ l e u + k v e tl ) 
+c 3 {k^qv + Kq^) + c^q v 
+C5 (q^ (35) 



where the c$'s are complex coefficients, and 

kft = (k , k), 
q^ = (k , -k), 

with 



k 2 =q 2 = 0, 

k-q-' 1 "^ 2 
k-e-- 



(ko) 
q ■ e 



= -l. 



(36a) 
(36b) 



(37a) 

(37b) 
(37c) 
(37d) 



The expansion (34)-(35) is the most general one 
for both vectors and symmetric rank-2 tensors and 
must be supplemented by the sources constraints 
(|30p . Accordingly, the positivity of the residue of 
the propagator is assured if 

3mRes(n| fe 2 =0 ) > 0. (38) 



III. TESTING THE EFFICACY AND 
QUICKNESS OF THE PRESCRIPTION 

In order to explicitly illustrate the generality and sim- 
plicity of the proposed method we analyze in the fol- 
lowing the unitarity of some higher-derivative electro- 
magnetic (gravitational) models enlarged by both Chern- 
Simons and higher order Chern-Simons terms. We also 
comment, in passing, on some interesting and remarkable 
properties of these systems. 



A. Higher-derivative electromagnetic models 

We begin our study by checking the unitarity of the 
Lee-Wick-Chern-Simons model which is defined by the 
Lagrangian, 



£lwcs = --F» V F^ - JLf^DF'*' + ^A^Ap, 

(39) 

where F^ v = d^A u — d^A vl and m (p) is a parameter 
with dimension of mass. 

Now, writing the Lagrangian above in the form (26) 
we promptly obtain the expression for the wave operator 
in momentum space, namely, 

<V = {^-k 2 + ^) V - «MW fcP - (40) 

With the help of the identities listed in the Appendix 
A, the wave operator (|40[) may be expanded in the 3D 
spin projection operators basis, as follows 

<v = E a ( J )^f( J U' ( 41 ) 
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where 



o(l) 



o(0) = 0, 



^h + ^Vk 2 

m " 





MA 



(42) 



(43) 



It is worth noticing that the spin-0 sector is completely 
degenerate, which is fully expected since the model ([31?)) 
has a gauge symmetry 



A' = A, 



8A, 



(44) 



The term 5A^ can be easily obtained by noticing that it 
can be associated with the right null eigenvalues of the 
matrices of the coefficients V- 



{R,n) 



i,J,0 



(45) 



This result applies to every independent value of j and 
n. Accordingly, we get 

SA, = d, {duf") , (46) 

where f u is an arbitrary function. 

Interestingly, the gauge symmetry of the model inhibits 
the propagation of the spinless mode; as a consequence 
of this symmetry, there appears a source constraint that 
prevents the propagation of this unphysical state. It is 
trivial to see that the general expression (|3"0|) reduces now 
to 



0. 



(47) 



which is nothing but the familiar source conservation re- 
lation. 

The inverse matrix of the spin-1 sector, on the other 
hand, reads 



where A 
(K-2K 



1 / -fc 2 



K - ma/P 
-fc 2 + 4 



fc 2 - 



n 2\ 2 k 2 



I* 




M 



k 2 is a quartic polynomial 

in k 2 . As a result, it has four roots: one massless pole, 
and three massive ones which we shall call m 1; m 2 , and 
TO3, respectively. Therefore, as far as the nature of the 
roots are concerned, there exist precisely four distinct 
cases to be dealt with (see Fig. 1): (1) fi — (Lee- Wick 
electrodynamics), (2) < /i 2 < (In this case there 

are three real positive masses.), (3) > ^S- (Here 
there are necessarily two complex roots, implying in the 
existence of tachyonic excitations.), and (4) p? = 
(In this case there appears a double pole; as a result, 
the model is non-unitary.). We discuss in the following 
only the physical models, i.e., cases (1) and (2). 




11=0, — < n < 

27 



2 4 m 
(1 > 

27 



Figure 1: Polynomial function A ^ - versus fc 2 , where A (k 
refers to the denominator of the propagator (|48H . 

1. /i = (Lee-Wick Electrodynamics) 
The matrix of the coefficients is now given by 



a(l) = 







-fc 2 + K 

while its inverse can be written as 



(49) 



a 1) 



1 



{k 2 -m 2 )k 2 \ 1 



1 



(50) 



Consequently, the absence of tachyons and ghosts in 
the model is subordinated, respectively, to the following 
conditions 



mf > 0, 



(-ljtr^l.m 2 )- 1 !^^ = -1, 



(51) 



(52) 



which clearly shows the presence of a non-tachyonic mas- 
sive ghost in the system. 

For the massless pole, the constraint fc M 5 M = allows 
us to write the saturated propagator as 



H ~ (k 2 ~m 2 )k 2tS *^- 



(53) 



Expanding now the current 5^ in the momentum basis 
h yields 



SmRes(n| fc 2 =0 ) = |ci| > 0, 



(54) 



s 



which allows us to conclude that the massless mode does 
not violate the unitarity. 

The wrong sign of Eq. (53) indicates an instability 
of the theory at the classical level. From the quantum 
point of view it means that the theory in non-unitary. 
Luckily, these difficulties can be circumvented. Indeed, 
the classical instability can be removed by imposing a fu- 
ture boundary condition in order to prevent exponential 
growth of certain modes. However, this procedure leads 
to causality violations in the theory (36|; fortunately, this 
acausality is suppressed below the scales associated with 
the Lee- Wick particles. On the other hand, Lee and Wick 
argued that despite the presence of the aforementioned 
degrees of freedom associated with a non-positive definite 
norm on the Hilbert space, the theory could nonetheless 
be unitary as long as the new Lee- Wick particles ob- 
tain decay widths. There is no general proof of unitarity 
at arbitrary loop order for the Lee- Wick electrodynam- 
ics; nevertheless, there is no known example of unitarity 
violation. Accordingly the Lee- Wick electrodynamics is 
finite. Therefore, we need not be afraid of the massive 
spin-1 ghost. We remark that recently a quantum bound 
on the Lee- Wick heavy particle was found that is of the 
order of the neutral vectorial boson found in nature 13 71 . 



implying that m\ < m 2 , m\ > m 2 , and to 2 < m 2 , which, 
of course, contradicts the original assumption that mi < 
vni < TO3. Thence, we come to the conclusion that this 
model is plagued by ghosts. 



Lee-Wick- Chern-Simons model enlarged by a higher 
derivative Chern-Simons extension 



Another interesting model can be built up from the 
Lee-Wich-Cern-Simons system by adding to the La- 
grangian (|39[) the higher derivative Chern-Simons exten- 
sion proposed by Deser and Jackiw [38] 



£ecs = -e^ p DA^A p . 



(58) 



2. < fj. < =p- 



Here A (fc 2 ) = (fc 2 



(fc 2 



k 2 /r 



where mi, 1712, and 7713 are the three real positive roots 
of A. We assume without any loss of generality that 
mi < m2 < ni3 (see Fig. [I]). On the other hand, the 
relations (— 1) tvA(l, m 2 )\ k 2 =m 2 > (i = 1,2,3,) lead to 
the following inequalities 



m 1 



m 2 



m 3 





(m\ 


— m 2 ) 




(m\ 


-m 2 2 
(ml 


)(mf- 
— m 2 ) 


m 2 ) 


(m\ 


— m\ 
(m§ 


)(ml- 
— m 2 ) 


mf) 


(m 2 


— m 2 


) ( m i - 


m\) 



<o, 



<o, 



< 0, 



(55) 
(56) 
(57) 



Let us then check the unitarity of this curious model. 
Starting from the wave operator in momentum space, 



O, 



w = (-fc 2 + ^) + Afc 2 ) W fcP , (59) 



it is straightforward to show that the spin-1 matrix of 
the coefficients and its inverse are respectively given by 



ail) 



-k 2 + ^ + ^ + Xk 2 )V¥ _\ 

-fc 2 + i4-( M + Afc 2 )^ ) ' 



(60) 



±r - (fx + Xk 2 ) Vk 2 

_fc 2 + ^ + ( M + Afc 2 )v / fc 2 

I 



(61) 



where 
A = 



An analysis similar to that done in Sec. Ill A tells us 

fc 2 . (62) 
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that the addition of the higher derivative Chern-Simons 
extension does not improve the non-unitarity of the Lee- 
Wick-Chern-Simons model. Neither does it cure the non- 
unitarity of the Lee- Wick system. 

B. Higher-derivative gravitational models 

Higher-derivative gravity augmented by a Chern- 
Simons term is defined by the Lagrangian 

£ = ^{aR + pR^R^ +lR 2 ) + %Ccs, (63) 
where 

c cs = s^r^ (a v r;, + lr: x r; K )j (64) 

is the Chern-Simons term and a, /3, 7, and /i are arbitrary 
coefficients. 




Due to the gauge symmetry of the model the spin-0 ma- 
trix above is evidently non-invertible which is translated 
into the usual source conservation constraint on the grav- 
itational sources k^S^ = 0. 



In the weak field approximation (g^ = rj pv + h pu ), this 
Lagrangian reduces to 

c m = ^--h^nh^ + lhnh-hdn&hr 

+h^d^d p h p ,) + j(h^U 2 h pu + hU 2 h 

-2h> iV Ud ll d p hJ> + 2h^d p d v d p d a h prJ ) 
+ 1 {hU 2 h - 2hnd fl d v h> MJ + h^d^dpd^h"^ 
+ tl h » £ ^p dx [p hvp _ dud a h°) . (65) 

We have now all the ingredients to compute the wave 
operator O pv ^ pa and expand it in the appropriate degree 
of freedom basis with the aid of the identities collected in 
the Appendix A. The resulting matrices of the coefficients 
concerning this expansion are 



:)fc2 o), ( 6 6) 

k 2 \ 

(a + pk 2 -n^jk 2 J ' 

I 

On the other hand, the inverse (^4 (O)^ 1 ) of the largest 
nondegenerate matrix extracted from a (0) , as well as the 
inverse of a(2)~ 1 , can be written as 



J 




((3/3 + 87) k 2 -a 




A(O)- 1 = 



[(3/3 + 87) k 2 -a)k 2 ' 



(68) 



a(2y 1 



(a + (3k 2 ) 2 - p?k 2 k 2 



a + /3k 2 - iiVk 2 






a + /3k 2 + nVk 2 



(69) 



Using the constraints (|52"|) on the matrices (|6"5)) -(|6"9" ]l . 
the following relations for the parameters are obtained 



Spin-2 : a < 0, /3 > 0; (70) 
Spin-0 : a > 0, 3/3 + 87 > 0. (71) 



Accordingly, for arbitrary values of the parameters the 
model is non-unitary. Nevertheless, there exists a way 
of circumventing this difficult: all we have to do is to 
prevent the propagation of the massive spin- mode by 
choosing 3/3 + 87 = 0. Remarkably, this is precisely the 
constraint utilized by Bergshoeff, Hohm, and Townsend 
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(BHT) in the construction of their model [l9| . Another 
alternative is to inhibit the propagation of the massive 
spin-2 mode by setting /3 — 0. As a consequence, 



A(oy 



(72) 



If we take into account that the absence of tachyons and 
ghosts requires respectively that — > and 7 > 0, we 

come to the conclusion that the model aR + ^R 2 is uni- 
tary if a and 7 are positive. 

For the massless poles, one must use the original ex- 
pression pip for the propagator in order to compute the 
residue. The constraints satisfied by the sources allow 
us to handle correctly the singularities. Using such con- 
straints and discarding terms that do not contribute to 
the residue, yields 



n 



1 



iS 



ak 2 



2 (VppVvcr 
X 



s pa . 



(73) 



Using a suitable basis for the expansion of the sources 
in momentum space, we arrive at the conclusion that this 
expression vanishes identically, which clearly shows that 
the massless mode is non-propagating. 

The preceding analysis seems to indicate the existence 
of two unitary higher-derivative gravity models in 3D: the 
BHT and the aR + ~fR 2 systems. Actually, only the BHT 
model can be really considered a higher-derivative gravity 
system. Indeed, this model contains fourth-derivatives of 
the metric, while the pure scalar curvature system is con- 
formally equivalent to Einstein gravity with a scalar field 
(3{|, which means that despite having fourth derivatives 
at the metric level the aR + 7-R 2 model is ultimately 
second-order in its scalar-tensor version. 



Therefore, up to now, "New Massive Gravity" [2fl, |21| 
is the only higher-derivative gravity model in 3D that 
is both perturbatively renormalizable and unitary in flat 
space [40j. Most likely the full model is non renormaliz- 
able since it improves only the spin-2 projections of the 
propagator but not the spin-0 projection. We point out 
that the description of gravitational phenomena via the 
BHT model does not lead to some really bizarre results as 
in the usual 3D gravity (lack of a gravity force in the non- 
relativistic limit, gravitational deflection independent of 
the impact parameter, complete absence of gravitational 
time dilation, no time delay) . Actually, in the framework 
of New Massive Gravity, short-range gravitational forces 
are exerted on slowly moving particles; besides, the light 
bending depends on the impact parameter, as it should 
fill ]. And more, both time delay and spectral shift do 
take place in the context of the alluded model (42j | . 
1. Higher-derivative-Chern-Simons gravity enlarged by the 
Ricci- Cotton tensor 



For gravity theories there is also the possibility of the 
construction of a higher derivative Chern-Simons exten- 
sion, the so-called Ricci-Cotton term, which is defined by 
the Lagrangian 



£ RC = 



'I? 



D„R p °. 



(74) 



We remark that models including this term were recently 
investigated by Bergshoeff, Hohm, and Townsend in their 
researches on higher derivatives in 3D gravity and higher- 
spin gauge theories [43]. 

Let us then probe the unitarity of higher-derivative- 
Chern-Simons gravity augmented via the Ricci-Cotton 
term. It is curious that this term only alter the spin-2 
sector of this model. The matrix of the coefficients and 
its inverse are now given by 



a (2) 



/3k 2 - (n + \k 2 ) 




13k 2 



Xk 2 



k 2 k 



(75) 



a (2)- = ! 



+ (3k 2 + (n + Xk 2 ) 




a + (3k 2 - (n + Xk^Vk 2 



(76) 



where 

A = [-A 2 fc 6 + ((3 2 - 2fiX) fc 4 + (2a(3 - ^ 2 ) k 2 + a 2 ] k 2 . 

A cursory glance at the equation above is sufficient 
to convince us that the model at hand can describe at 
most three massive particles. Proceeding in the same 



way as we have done in Sec. Ill A we conclude that the 
addition of the Ricci-Cotton term is not a good therapy 
for curing the nonunitarity of higher-derivative-Chern- 
Simons gravity. 
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IV. CONCLUDING REMARKS AND 
COMMENTS 

We have devised an easy procedure to check the unitar- 
ity of 3D models based on a new class of spin projection 
operators. The great importance of these operators re- 
sides precisely in the fact that they form a complete set 
for Lorentz preserving (PV and PC) models. In other 
words, they obey a completeness relation. There are, 
however, some discussions about the validity of using an 
"incomplete" set of operators, i.e, a set that does not have 
all the elements it needs, to expand the propagator. In 
general, it is alleged that the operation of multiplication 
of the operators of the alluded incomplete set obeys the 
axiom of closure. The fact that a set is closed under a 
given operation does not imply of course that it is com- 
plete. That is, the term closure cannot be used as a 
synonym for completeness; indeed, any complete set is 
necessarily closed but the reciprocal statement is false. 
A simple example is sufficient to clarify this point. Con- 
sider, in this vein, three-dimensional PC gravity mod- 
els. In this case, as we have already commented, the 
appropriate basis for computing the propagator is that 
whose elements are the well-known Barnes-Rivers oper- 
ators. Expanding the wave operator in this basis, we 
obtain in momentum space 



O = x 1 P(l)+x 2 P(2)+x s P(0 s )+x LJ P{0 LJ ) 

+^P(0 S ")+^ S P(0" S ). (77) 
Consequently, the corresponding propagator is given 



by 



o- 1 = 



— P(2) + — P(l) + - 

x u P(O s ) + x s P(0") - x su P(O s ") 
-x us P(0" s ) 



(78) 



Now, if S = {P(1),P(2),... ,P((T S )}, then S' = 
{P(l), P(2), P(0 S ), P(0")} is a subset of S which is closed 
under the same operation of multiplication as that con- 
cerning S; in addition, the elements of 5" obey the rela- 
tion 



P(l) + P(2) + P(0 s ) + P((n = 5, 



(79) 



which is nothing but the decomposition of unity. The 
point, nonetheless, is that it is sometimes claimed that 
(79) is a completeness relation for the operators at hand. 
If this were so, we would arrive at the wrong conclusion 
that 5" should necessarily be a basis for performing our 
computations; as a consequence, the propagator for the 
PC gravity models would assume the form 



1 



1 



1 



1 



O-' = -P(l) + —P(2) + -P(0 S ) + — P(O-). 



X2 



(80) 



Comparing (78) and (80) we come to the conclusion 
that these expressions coincide only and if only x su> — 
Xujs — 0. Nevertheless, these coefficients cannot be zero. 
Indeed, since the PC gravity models are gauge invariant, 
we have to add to the Lagrangian of the model a gauge- 
fixing Lagrangian (£ g f) so that the resulting wave oper- 
ator can be inverted. Choosing for this purpose, without 
any loss of generality, the de Donder gauge and taking 
into account that its linearized version can be written as 
follows 



(81) 



where 7^" = d v h^ v — ^d^h, we promptly obtain in mo- 
mentum space 



O gf (k) = 



2 



-P(l) + -P(0 S ) + -P(0 W ) 



^P(0^ 
4 



(82) 



which clearly shows that both x SUJ and x us are differ- 
ent from zero. In other words, the operator wrong = 
ZiP(l) + x 2 P(2) + x s P(0 s ) + ir w P(0 w ) is obviously no 
invertible. Suppose, however, that one stubbornly in- 
sists that both expressions for the propagator are cor- 
rect due to the fact that for physical problems in which 
the propagator (78) is contracted with conserved exter- 
nal currents (JO~ 1 J,kJ = 0), both operators P(0 S ") 
and P(0 WS ) do not contribute for the final result of the 
calculations. Again, it is trivial to show that this argu- 
ment is fallacious. In fact, a straightforward computation 
leads to the following results 



JO- 1 J = J 



JO' 1 T — J 



1 



P(2) 



X2 

— P(2) + — P(0 S ) 

X 2 



X s X(j} Xg^X^s 
J. 



-P(0 S ) 



In summary, only operators that form a complete set 
can be used to attaining the propagator. 

Another point that deserves to be discussed is whether 
or not the nonunitary disease that affects some three- 
dimensional models could be cured by the addition of 
a Chern-Simons term to the system; of course, we are 
not excluding from our considerations the possibility 
of enlarging the model via a higher derivative Chern- 
Simons extension or even through the simultaneous ad- 
dition of Chern-Simons and higher derivative Chern- 
Simons terms. Before answering this important question, 
let us briefly comment, in parenthesis, on the origin of 
the widespread habit of augmenting the Lagrangian of a 
given model through a Chern-Simons term. For the sake 
of brevity, we restrict our analysis to three-dimensional 
gravitational models. Everything started when it was 
found out that the solution to the triviality problem of 



general relativity in (2+l)D could be cured by simply 
adding a topological Chern-Simons term to the system. 
The resulting model describes a non-trivial gravity theory 
with a propagating, massive, spin-2 mode [7|. Later on it 
was considered another way out of the triviality problem 
of 3D gravity: the addition of higher-derivative terms 
to the system [44] ; unfortunately the resulting models 
are nonunitary [45] . On the other hand, it was claimed 
that the addition of a Chern-Simons term to the previous 
model would cure its non-unitarity (46| . This was proved 
afterwards to be incorrect |47j. After this digression, let 
us respond the question we have raised above. As we have 
seen in Section III, nonunitary higher- derivative elec- 
tromagnetic (gravitational) models do not become uni- 
tary systems by simply augmenting them through Chern- 
Simons terms. Neither do they become unitary by enlarg- 
ing them via a higher derivative Chern-Simons extension. 
It is amazing, nonetheless, that are some examples in the 
literature of unitary systems whose unitarity is spoiled by 
the addition of Chern-Simons terms [48l. |49|. Therefore, 
in some cases, the coexistence between the topological 
term and higher-derivative theories is conflicting. 



satisfied by them. 



1. Vector field operators: A — A 



a. Spin-0 Sector 



b. Spin-1 Sector 



p£ a (+iL„ = /v 



c. Identities Among the Operators 
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Appendix A: Projection Operators and Tensor 
Relations 



In this appendix, we gather the degree-of-freedom op- 
erators constructed in Sec. Ill Al and some useful identities 



• P AA = <V = Pii + Ptf (- 

d. Tensorial Identities 

= P AA (0)^ + P AA (l)^ 
k p k v =k 2 P AA (0)^ 



J 



2. Rank-2 Symmetric Field Operators: h — h 



e. Spin-0 sector 

• P/f {0 s ) p „. piJ = \0 pv d pa , 

• P22 {® U )uv;p<T = UfivWpa, 



phh (QSU\ 
P^ (0" 



' pv.pa 

f. Spin-1 Sector 

• Pit 



livypcy 



^^livOpcr- 



\ {Pup^va + Pvp^po + PfujLdup + p va U)pp) , 
\ {o~ P p^ua + 0~vp^^a + 0~ prJ UJ up + 0~ l/a LO pp ) , 



Pl2 (± l )^pa = \ E -rr,K (p p <J V p^ a + pl^UJ^ + p^O^vp + pl^UJ w ) , 

P21 (Tl)^ ;pCT = \e T7]K {olp T p u av + olp T a u pv + a^p^w^ + a^p T a uj pp ) 
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g. Spin-2 Sector 

• Pii (+2) M1/;/M7 = Piippva, 

ft. Identities Among the Operators 

• p Wi W„ - \ ( + V«w + + - Pu (+i) + ^ (-i) 

• p Wi (2)„ = I + e^e vp e^e pa ) = (+2) + (-2) 

i. Tensorial Identities 

S^pa - £ i^pVua + VnaVvp) = (2) + (1) + A? (0) + ^ (0) 

= 2P?i (0 s ) + V2P& (0) + (0) + (0) 

*V.Mp* + fcpfcaV = ^ fc2 (P12 1 (0) + P2? (0)) + 2fc 2 P 2 fc 2 h (0) 
kpkpT\ V(7 + k^kfjfjjyp + hj/hpTjfxcf + hi/hcpfjfxp = 2k P (1) + 4& _P 2 2 (0) 
kpki/kpkfj — k P22 (0) 

2VF (P/f (+2) - f* fc (-2) - (±1) + P 2 'f ( T l)) 

{^K,pnkiik a + s K p U k^k a + s Ka ^k u k p + s Kav kpkp) k — 2k Vfc 2 ( P 12 (il) "I" P21 (~Fl)) 
Here the /zz/; per indices of the operators P^ {J) pv . pa were omitted. 
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